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ABSTRACT 
Many ecological studies and conservation 
management plans employ noninvasive scat 


sampling based on the assumption that species’ scats 
can be correctly identified in the field. However, in 
habitats with sympatric similarly sized carnivores, 
misidentification of scats is frequent and can lead to 
bias in research results. To address the scat 
identification dilemma, molecular — scatology 
techniques have been developed to extract DNA from 
the donor cells present on the outer lining of the scat 
samples. A total of 100 samples were collected in the 
winter of 2009 and 2011 in Taxkorgan region of 
Xinjiang, China. DNA was extracted successfully from 
88% of samples and genetic species identification 
showed that more than half the scats identified in the 
field as snow leopard (Panthera uncia) actually 
belonged to fox (Vulpes vulpes). Correlation between 
scat characteristics and species were investigated, 
showing that diameter and dry weight of the scat were 
significantly different between the species. However it 
was not possible to define a precise range of values 
for each species because of extensive overlap 
between the morphological values. This preliminary 
study confirms that identification of snow leopard 
feces in the field is misleading. Research that relies 
upon scat samples to assess distribution or diet of the 
snow leopard should therefore employ molecular 
scatology techniques. These methods are financially 
accessible and employ relatively simple laboratory 
procedures that can give an indisputable response to 
species identification from scats. 


Keywords: DNA analysis; Snow leopard; Scats; 
Noninvasive genetics; Carnivore 


INTRODUCTION 


Information gathered from scats, such as diet, distribution, 
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abundance and community dynamics, is widely used in many 
ecological studies and conservation management plans and 
requires reliable identification of scats (Gibbs, 2000; Long et 
al, 2008). However the identification of species from field 
signs alone is not always accurate and misidentification of 
scats is frequent, in particular in habitats with sympatric, 
similarly sized carnivore species (Farrell et al, 2000). 
Traditional methods of identifying scats in the field include size 
(length and/or diameter), shape (segmented, pointed ends, 
etc.), color (black, brown, white, etc.), pH or smell, in addition to 
the co-occurrence of other signs at the site such as pugmarks, 
scratches or hair (Danner & Dodd, 1982; Green & Flinders, 
1981; Jackson & Hunter, 1996). However, these field methods 
have proven unreliable for the following reasons: body size can 
vary greatly within species which affects scat dimensions (e.g. 
between juvenile and adult, and male and female), and other 
species may investigate the scat sites of the target species, 
confusing identification by leaving their signs. i 

To address the scat identification dilemma, molecular 
scatology techniques have been developed to extract DNA 
from the donor cells present on the outer lining of the scat 
samples ( Foran et al, 1997; Höss et al, 1992). Target DNA, 
that enables species identification, is then amplified using 
polymerase chain reaction (PCR) based methods. These 
techniques have been employed for scat recognition in 
numerous studies and for a variety of species, for example from 
identification of European brown bears (Ursus arctos) (Kohn et 
al, 1995), to sympatric confamilial groups (Mills et al, 2001), and 
discerning species in an entire community of carnivores 
(Fernandes et al, 2008) or even telling apart hybrids (Adams et 
al, 2003).The use of molecular scatology has also revealed and 
confirmed that species identification of scats in the field is often 
inaccurate (Davison et al, 2002; Farrell et al, 2000; Harrington 
et al, 2009; Prugh & Ritland, 2005; Zuercher et al, 2003). 
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Snow leopards are very elusive, extremely rare and inhabit 
harsh environments (Jackson & Ahlborn, 1988; Sunquist & 
Sunquist, 2002) making research through direct observation 
almost impossible. Research on this species therefore relies 
on gathering information in the absence of direct contact, 
primarily through non-invasive techniques (Ma et al, 2011; 
Wei et al, 2001; Xu et al, 2005). One of the main survey 


techniques used is scat collection with the following objectives: 


to determine prey preference (Bagchi & Mishra, 2006; Liu et 
al, 2003; Oli, 1993; Shehzad et al, 2012; Wang et al, 2014), 
to estimate population size ( Ale et al, 2007; Fox et al, 1991; 
Hussain, 2003; Karmacharya et al, 2011; McCarthy et al, 


2008), to analyse community structure (Lovari et al, 2013), etc. 


In most of their ranges, snow leopards co-occur with other 
similarly sized carnivores, such as wolves, foxes and jackals 
(Canis aureus), potentially creating uncertainty in the field 
identification of scat. It is not surprising that evidence of extensive 
scat misidentification confirmed by molecular techniques has 
been reported for snow leopards in previous publications. 
Janečka et al (2008) first observed a high level of scat 
misidentification, ranging from 35% in Ladakh, India to 54% in 
South Gobi, Mongolia. Of 71 samples collected in the study by 
Karmacharya et al (2011) in Nepal, 42% were misidentified as 
snow leopard in the field, but in fact belonged to other carnivores. 
Anwar et al (2011) in Pakistan found 52% of scats to be correctly 
identified as snow leopard while the rest were from other 
sympatric species. Likewise Shehzad et al (2012) in Mongolia 
correctly identified in the field 43% of putative snow leopard scats 
as confirmed by DNA analysis, while 57% were excluded from 
the study because they actually belonged to other species. 
Despite misidentification being a source of significant bias, 
scat sampling still remains a fundamental tool in the study and 
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conservation of snow leopards, so it is important to address 
and reduce the described inaccuracies. 

This study aims to: (1) assess the extent of misidentification 
in the field for snow leopard scat; (2) determine if variables 
associated with the scat, such as morphological characteristics 
of the feces, are correlated with species identification; and (3) 
understand if accurate field identification of snow leopard scats 
is possible, if so establish more specific field collection protocols 
that include morphological features statistically relevant for 
successful species identification of scats. 


MATERIALS AND METHODS 


Study area and sample collection 

Scats were collected from February 22, to March 12, 2009, and 
March 12 to April 12 in 2011 in Taxkorgan Nature Reserve 
(TNR; E74°30’-77°00’, N36°38’-37°30’), located in the east 
plateau of the Pamir Mountains, Xinjiang Uygur Autonomous 
Region, China. TNR has a mean elevation of approximately 
4000 m and it is characterized by a cold desert climate, with 
long very cold winters. The average temperature during the 
survey months is 1.03 °C, with 3.27 mm average precipitation 
per month. The carnivore guild includes snow leopards, wolf 
(Canis lupus), red fox (Vulpes vulpes) and lynx (Lynx lynx). 
Pallas' cat (Otocolobus manul) may also be present, but 
extended surveys using camera traps have never recorded in 
our study sites. 

Data for this study were collected in two sites, Mariang and 
Mazar; both sites had minimum altitudes above 3 000 m. The 
study sites straddled the southeastern and eastern 
boundaries of the TNR respectively, with some transects 
extending as far as 40 km outside the reserve (Figure 1). In 
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Figure 1 Location of Taxkorgan Nature Reserve including the two study sites (Mariang and Mazar) 
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2009, 24 transects with a total length of 242.7 km were di- 
vided between the two study areas, with 13 transects covering 
204.1 km in Mazar, 11 transects totalling 38.6 km in Mari- 
ang. In 2011, 40 transects with a total length of 131.6 km were 
conducted in Mariang. Each transect was located within a 
distinct valley or along a ridgeline. In the field, scat samples 
were attributed to a certain species according to 
morphological characteristics (Table 1), the presence of other 


discerning signs (e.g. pugmarks, scratches, urine etc.) and 
the scat’s location (e.g. under a cliff, open area, on trail etc.). 
These parameters were used in combination as guidance to 
identify the correct donor species. This was left to some 
extent to the collectors’ discretion and decisions could be 
based on experience and individual opinion. The field team 
was composed of academic professionals specialized in 
carnivore studies and local Nature Reserve staff. 





Table 1 Scat morphological parameters used in the field as guidance to identify the species. 

Species Color Shape Number of segments (n) Length (cm) Diameter (cm) 
Snow Leopard Dark Long pointed tail <3 5-6 >2 

Lynx Dark Long pointed tail <3 4-5 ~2 

Wolf/dog White or grey Round end, no tail <3 5 >2 

Red Fox Dark, grey old samples Round end, short pointed tail 2-5 2-4 <2 


Each sample at time of collection was placed individually in 
a labeled zip lock plastic bag to avoid contamination. To 
prevent degradation the samples were air-dried in the field by 
placing the open bags in a ventilated, cool and dry 
environment away from direct sunlight. In the laboratory, 
samples were then transferred to a freezer at -4 °C and 
finally to -80 °C for long-term storage. The morphological scat 
data recorded included number of segments (S), mean 
segment length (ML, cm), total length of scat (TL, cm), mean 
diameter (MD, cm) and weight after drying (DW, g). Original 
weight at collection was also noted but not used in these 
analyses, as it is greatly dependent on the age of the scat 
(fresh scats are heavier due to a higher water content, while 
old scats are lighter) and location of the scat (exposure to ice 
and snow may increase water content of old scats). Weight 
after drying was preferred as it does not present this 
constraint and was considered more appropriate. 

Scat samples (number of samples in brackets next to the 
species Latin name in bold) were assigned to reference 
species based on nucleotide diversity <0.03 and node 
bootstrap value >90%. 


DNA extraction and species identification 

Laboratory analysis was performed at the Key Laboratory of 
Animal Ecology and Conservation Biology (Institute of 
Zoology, Chinese Academy of Sciences, Beijing). The DNA 
extraction was performed using the QlAamp DNA stool minikit 
(Qiagen). A 146 bp region of the mitochondrial cytochrome b 
gene was amplified by PCR using carnivore specific primers 
from Farrell et al (2000) (5'-AAACTGCAGCCCCTCAGAATG 
ATATTTGTCCTCA-3' and 5'-TATTCTTTATCTGCCTATACAT 
RC ACG-3'). It is not possible to distinguish between wolf and 
domestic dog using this primer due to their close phylogenetic 
relationship, so they will be referred to as “wolf/dog” in the rest 
of this paper. Amplifications were conducted following the 
protocol developed by Janecka et al (2008), although changes 
to the volume of reagents (total 50 uL instead of 10 uL) and 
thermo cycling conditions were made as follows: 5x 
PrimeSTAR buffer 10 uL, 4 uL of dNTP mixture (2.5 mmol/L), 
forward primer 1 uL (10 pmol/L), reverse primer 1 uL (10 
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umol/L), PrimeSTAR HS DNA polymerase 0.5 uL (2.5 U/uL), 
DNA extract 4 uL and DNA-free water 29.5 uL. The PCR 
conditions included an initial denaturing step of 94 °C for 1 
min, followed by 40 cycles of 94 °C for 30 sec, 54 °C for 30 
sec, 72 °C for 30 sec, and a final extension step of 72 °C for 2 
min. PCR products (3 uL amplified DNA plus 1.5 uL of 
loading dye) were fractionated on a 2% agarose gel, stained 
with ethidium bromide, and visualized under ultraviolet light. 
Primers were provided by Invitrogen (Beijing, China) and all 
reagents by TaKaRa Biotechnology Co. (Dalian, China). PCR 
products were sent for sequencing to SinoGenoMax Company 
Ltd. (Beijing, China). Sequences obtained from one strand 
were examined in 4PEAKS version 1.7.2 (©2006 Mek& 
Tosj.com) and submitted to a BLAST search (Madden, 2002) 
to be compared with entries in GenBank sequences in order 
to identify the species for each sample (Zhang et al, 2000). 
Sequence matches with an E-value equal to 0.0 and 
maximum identity value 290%, were considered positive 
species identifications (DeMatteo et al, 2014; Keehner, 2009; 
Naidu et al, 2011; Rozhnov et al, 2011). In addition, a 
phylogenetic tree was constructed with reference species 
(downloaded from GenBank with accession numbers: 
KJ637144, JQ003577, AB303951, AB194817, EF689046, 
AY928671, KF990330, EF551002, AF053050, JF357970, 
JF357968, EF551004, KF661088) to double-check species 
identification, using the neighbor-joining algorithm based on 
the Kimura 2-parameter model (Kimura, 1980) with MEGA 
v6.06 (Tamura et al. 2013). Node support was evaluated 
using 1 000 bootstrap replicates. 


Variables associated with species identification 
The genetic species identification results were used to 
calculate field identification error for snow leopard feces (i.e. 
the percentage of scats that were misidentified as snow 
leopard in the field). A confusion matrix (or contingency table) 
was built to summarize the results of the misidentification for 
all collected samples belonging to the carnivore guild. 

A one-way ANOVA test was used to look for correlation 
between the presumed snow leopard samples (categorical 
independent variables) and the morphological variables 


(continuous dependent variables) associated with the scat as 
previously done by Anwar et al (2011). For those variables 
that were found to be significant, a post-hoc test of least 
significant difference (LSD) to explore all possible pair-wise 
comparisons of means comprising a factor using the 
equivalent of multiple t-tests. 

To analyze the relationship between the variables 
associated with the presumed snow leopard scat and the 
actual species, principal component analysis (PCA) was used. 
This method is a true eigenvector-based multivariate analysis, 
which reduces the effective dimensionality of a multivariate 
data set by producing linear combinations of the original 
variables that summarize the predominant patterns in the data 
(Peres-Neto et al, 2003). 

All analyses were performed using software R 2.11.1 (R 
Development Core Team, 2013). The total number of samples 
included in the analysis varied as some samples had missing 
values (i.e. in 2009 total length, mean segment length and 
number of segments were not recorded) and could not be 
used in certain statistical tests. 


RESULTS 


Identification of species from scat 

During the survey 100 scats were collected and were 
attributed to a species in the field: 51 presumed snow leopard 
scats, 40 presumed fox scats and 9 could not be identified 
(Figure 2). No wolf/dog samples or other carnivores were 
categorized in the field. DNA was successfully extracted, 
amplified and sequenced for species identification from 88% of 
scats and all species were categorized with maximum identity 
value290% (Appendix 1, available online). Nucleotide diversity 
(number of nucleotide base substitutions per site averaged over 
all sequence pairs within each species) was 0.005 for snow 
leopard, 0.009 for fox and 0.024 for wolf/dog. Samples were 
assigned to reference species based on pairwise distance 
<0.03 and node bootstrap value >90% (Figure 3; Appendix 2, 
available online). Unidentified scats (12 samples) could either 
be from species for which the primers were not appropriate 
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Figure 2 Species identification of scat samples (“Field species 
ID”) in the field compared to results of genetic identification 
(“Genetic species ID”) 
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Figure 3 Condensed neighbor-joining tree with cut-off value of 
90% built using Kimura 2-parameter model 


(non-carnivore species) or represent samples with low DNA 
quality not suitable for genetic analysis. The discrepancy 
between species identification in the field and identification by 
genetic analysis is shown in Figure 2. 

Among the samples identified in the field as snow leopard 
scats (51 samples), the majority was genetically identified as 
actually belonging to fox (28 samples, 55%), 4 were wolf/dog 
(8%), while 16 were snow leopard (31%). The remaining 3 
were unidentified samples (6%). Therefore snow leopard scat 
identification error, calculated based on genetically verified 
scats for snow leopard scats in our study is 67%. 

Genetic analysis showed that scats categorized as fox in 
the field were correctly identified in all cases, except when 
DNA analysis was unsuccessful (Table 2, total 40 samples, 32 
confirmed fox and 8 were unidentified). A confusion matrix 
was built (Table 3), showing that while Type 1 errors were 
made in 32 occasions (scats identified as snow leopard in the 
field were not snow leopard scats), Type 2 errors didn’t occur 
(snow leopard scats were not present amongst scats 
identified in the field as other species). 

Table 2 Morphological characteristics (diameter and dry weight) 
for snow leopard, fox and wolf/dog. 

Species Total samples Diameter (cm) Dry weight (g) 
2.5+0.5 16.2415.1 
1.8+0.4 6.8+7.0 


2.440.6 14.848.3 


Snow leopard 16 
Fox 67 
Wolf/dog 5 


Table 3. Confusion matrix summarizing the results of 2009 and 
2011 field seasons 


Predicted class 








Snow Fox Wott! dog Unidentified Total 
leopard 
i Snow leopard 16 0 0 0 16 
& Fox 28 32 0 7 67 
= Wolfidog 4 0 0 1 5 
2 Unidentified 3 8 0 1 12 
< Total 51 40 0 9 100 


Type 1 errors were made (n=32) but no Type 2 errors were made 
when trying to identify snow leopard scats in the field. 
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Table 4 Results of snow leopard scat identification studies 





Author Country Total samples Snow leopard Fox Corsac fox Lynx Wolf/ dog Unknown ID error 
Janecka India 32 53% 19% NP 0 6% 22% 35% 
Janecka Mongolia 27 41% 48% NP 0 0 11% 54% 
Anwar Pakistan 95 52% 21% 3% NP 11% 13% 40% 
Karmacharya Nepal 71 27% 42% 31% 61% 
Shahzad Mongolia 203 43% Not analyzed 57% 
Present study China 51 31% 55% NP 0 8% 6% 67% 


Describe percentage of samples for each species and misidentification; ID error: percentage of feces that were incorrectly identified as snow 


leopard in the field; NP: not present in study site. 

Scat morphological variables of feces and_ their 
association with species 

The length of snow leopard scat samples was on average 
11.0+4.3 cm, ranging from shortest 3.6 cm to longest 18.1 cm. 
The average diameter was 2.1+0.5 cm and the weight after 
drying 16.2+15.1 g. The number of segments was between 1 
and 2 (2.2+1.2) with each segment having a mean length of 
5.9+3.1 cm. However the values for wolf/dog scats should 
only be used as reference value as the low sample size does 
not allow for statistical accurate results. 

One-way ANOVA tests showed that the diameter of scat 
was significantly different between snow leopard, fox and 
wolf/dog (F2,85)=18.13 and P<0.01), with significant difference 
between snow leopard and fox scats (LSD test, P<0.01) and 
wolf and fox (LSD test, P<0.01), but snow leopard and 
wolf/dog scats were not significantly different. In particular fox 
scats had a considerably smaller diameter (1.8+0.4 cm), while 
snow leopard and wolf/dog scats were larger and similar in 
size (2.5+0.5 cm and 2.4+0.6 cm for snow leopard and 
wolf/dog respectively) (Table 2). 

Dry weight also differed significantly between species 
(Fi2,85)=7.99 and P<0.01). Snow leopard and fox dry weights 
were significantly different (LSD test, P<0.01) but highly 
variable (16.2+15.1 g and 6.8+7.0 g for snow leopard and fox 
respectively), while wolf/dog scats were more consistent in 
weight (14.8+8.3 g) (Table 2). None of the other morphological 
variables tested were significantly different. The first two 
dimensions from the principal component analysis (PC1 and 
PC2) represented a cumulative proportion of variance of 79% 
(43% and 36% respectively). Snow leopard scat samples 
were distributed throughout the component space and no 
discriminating patterns could be identified. 


DISCUSSION 


The high misidentification rate of snow leopard scats indicates 
that field identification is problematic, so erroneous 
conclusions on the overlapping diet of these two species and 
incorrect prey preference were likely to be made if our study 
had proceeded without first employing molecular analysis. 
However, it is reassuring to know that amongst the feces field- 
identified as fox, no snow leopard samples were found, hence 
valuable samples were not lost. Other researchers studying 
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snow leopards have encountered similar issues with scat 
identification and their results are compared in Table 4. 
Overall, only a few of the total samples collected in these 
studies actually belong to the target species (Anwar et al, 
2011; Janecka et al; 2008, Janecka et al, 2011; Karmacharya 
et al, 2011; Shehzad et al, 2012), raising concern for research 
conclusions inferred from studies that did not employ 
molecular analysis. 

Identification error in our study decreased significantly 
between 2009 and 2011 (Pearson’s Chi-squared test with 
Yates’ continuity correction, x’=0.08, do=1, P=0.77), from 
71% to 65% respectively, suggesting that the collector’s 
experience may have a role in reducing error. This can be 
verified in the future subsequent to several sampling seasons 
being carried out. 

Very few samples were unsuccessfully genetically analyzed 
with only 6% of presumed snow leopard scats unidentified 
(and 12% of total scats). This is likely due to the dry climate 
and the constant freezing temperatures of the study site that 
helped to preserve high quality DNA. For this study simple air- 
drying was sufficient for mtDNA analysis and species 
identification. Other studies presented in Table 4 show higher 
percentages of unidentified samples. This is likely due to DNA 
degradation that could result from differences in collection 
method such as different season, sun exposure, storage 
method etc. (Stenglein et al, 2010). Temperature, humidity 
and microorganisms may affect final DNA quality and have 
proven to be detrimental for molecular analysis of scats, so it 
is advisable to seek the best protocol available to prevent 
degradation and contamination. Studies have highlighted that 
silica and ethanol are two effective ways of preserving 
samples (Conradi, 2006; Santini et al, 2007). 

Amongst the morphological variables tested, none could 
confidently be used as guidance to correctly discriminate snow 
leopard from fox or wolf scats. In common with previous results 
by Anwar et al (2011), statistical significance was found in scat 
diameter and dry weight. It was not possible, however, to define 
precise ranges for each species because of extensive overlap 
between the scat morphological values. Other variables related 
to the surrounding environment of the scat, such as presence of 
other signs, location and substrate type, may also play a role in 
correct species identification, and should be recorded and 
verified in future studies (Ma et al, 2005). 


Davison et al (2002) have emphasized that scat 
morphological identification methods need to be more 
rigorous when used in surveys and suggest using multi- 
evidence approaches involving a variety of methods to 
correctly identify species presence. This is also the case for 
snow leopard surveys as scat misidentification can lead to 
confusion about the species’ presence and risk of 
overestimating populations. It is important that research not 
rely solely on the identification of scats in the field, but scat 
identification be improved by including other techniques such 
as molecular analysis, scat detection dogs and camera 
trapping (Janecka et al, 2011; Long et al, 2007a; Long et al, 
2007b;). An extensive sampling survey required to collect 
sufficient snow leopard scat samples can incur high financial 
and time costs, so an additional cost of 70¥ (about US$12) 
per sample for molecular analysis is a reasonable price to 
make these surveys more accurate. These costs will likely 
decrease in the future and more cost-effective analyses can 
be used such as real-time PCR or PCR-RFLP (Cossios & 
Angers, 2006; Harrington et al, 2009; Mukherjee et al, 2010; 
Rodgers & Janečka, 2012). 

This preliminary study confirms that there is a high rate of 
misidentification of snow leopard scats in the field and that 
morphological characteristics of scats can’t be used to reliably 
differentiate between sympatric carnivore species. Therefore 
any research project that requires species identification (diet 
studies, sign surveys etc.) is advised to employ noninvasive 
DNA testing of scats in order to avoid serious bias in results. 
Genetic methods involve rather straightforward laboratory 
procedures, are relatively inexpensive and provide 
indisputable species identification of scats. 

It is also important to note that only a small portion of 
samples collected actually belonged to the target species, and 
since small sample size is already considered to be a problem 
when studying elusive animals such as the snow leopard, it is 
important to remind researchers that they must take this into 
further account and plan accordingly by increasing the study 
area, extending duration of surveys or using multiple non- 
invasive sampling methods. 

Finally, molecular species identification is not only a 
verification tool to be used in the laboratory. The results, in 
combination with other information on scat morphology, can 
also provide valuable feedback to field workers to improve 
collection guidelines and ultimately create effective 
conservation action plans. 


ACKNOWLEDGEMENTS 


This project was supported by State Forestry Administration of China (SFA), 
Snow Leopard Network, the Snow Leopard Trust, Panthera, Beijing 
Forestry University and the University of Oxford. Many thanks go to Dr 
Ming LI from the Institute of Zoology, Chinese Academy of Science, for 
access to a high quality laboratory and Dr Hui-Juan PAN from Beijing 
Forestry University for technical advice. We appreciate the great support 
from the Forestry Department of Xinjiang Uygur Autonomous Region, 
Wildlife Conservation Division of Kashi Region, Taxkorgan Nature Reserve 
Administration, Livestock Husbandry Bureau of Taxkorgan County, and 


government of Mariang village. We greatly thank our local guide Kemeti 
Taklashur for his tremendous contribution in the field. 


REFERENCES 


Adams JR, Kelly BT, Waits LP. 2003. Using faecal DNA sampling and GIS 
to monitor hybridization between red wolves (Canis rufus) and coyotes 
(Canis latrans). Molecular Ecology, 12(8): 2175-2186. 


Ale SB, Yonzon P, Thapa K. 2007. Recovery of snow leopard Uncia uncia 
in Sagarmatha (Mount Everest) National Park, Nepal. Oryx, 41(1): 89-92. 


Anwar MB, Jackson R, Nadeem MS, Janečka JE, Hussain S, Beg MA, 
Muhammad G, Qayyum M. 2011. Food habits of the snow leopard Panthera 
uncia (Schreber, 1775) in Baltistan, Northern Pakistan. European Journal of 
Wildlife Research, 57(5): 1077-1083. 

Bagchi S, Mishra C. 2006. Living with large carnivores: predation on 
livestock by the snow leopard (Uncia uncia). Journal of Zoology, 268(3): 
217-224. 


Conradi M. 2006. Non-invasive sampling of snow leopards (Uncia uncia) in 
Phu valley, Nepal. Master thesis. Available at http: //urn.nb.no/URN: NBN: 
no-19945. 


Cossios D, Angers B. 2006. Identification of andean felid feces using PCR- 
RFLP. Mastozool Neotrop, 13(2): 239-244. 


Danner D, Dodd N. 1982. Comparison of coyote and gray fox scat 
diameters. Journal of Wildlife Management, 46(1): 240-241. 


Davison A, Birks JDS, Brookes RC, Braithwaite TC, Messenger JE. 2002. On 
the origin of faeces: morphological versus molecular methods for surveying 
rare carnivores from their scats. Journal of Zoology, 257(2): 141-143. 


DeMatteo KE, Rinas MA, Argiielles CF, Zurano JP, Selleski N, Bitetti MSD, 
Eggert LS. 2014. Noninvasive techniques provide novel insights for the 
elusive bush dog (Speothos venaticus). Wildlife Society Bulletin, doi: 
10.1002/wsb.474. 


Farrell LE, Roman J, Sunquist ME. 2000. Dietary separation of sympatric 
carnivores identified by molecular analysis of scats. Molecular Ecology, 
9(10): 1583-1590. 

Fernandes CA, Ginja C, Pereira |, Tenreiro R, Bruford MW, Santos-Reis M. 
2008. Species-specific mitochondrial DNA markers for identification of non- 
invasive samples from sympatric carnivores in the Iberian Peninsula. 
Conservation Genetics, 9(3): 681-690. 

Foran DR, Crooks KR, Minta SC. 1997. Species identification from scat: an 
unambiguous genetic method. Wildlife Society Bulletin, 25(4): 835-839. 
Fox JL, Sinha SP, Chundawat RS, Das PK. 1991. Status of the snow 
leopard Panthera uncia in Northwest India. Biological Conservation, 55(3): 
283-298. 

Gibbs JP. 2000. Monitoring populations. /n: Boitani L, Fuller TK. Research 
Techniques in Animal Ecology, Controversies and Consequences. New 
York: Columbia University Press, 422. 

Green J, Flinders J. 1981. Diameter and pH comparisons of coyote and red 
fox scats. Journal of Wildlife Management, 45(3): 765-767. 

Harrington LA, Harrington AL, Hughes J, Stirling D, Macdonald DW. 2009. 
The accuracy of scat identification in distribution surveys: American mink, 
Neovison vison, in the northern highlands of Scotland. European Journal of 
Wildlife Research, 56(3): 377-384. 

Höss M, Kohn M, Pääbo S. 1992. Excrement analysis by PCR. Nature, 
359(6392): 199. 


Zoological Research 36(2): 72-78, 2015 77 


Hussain S. 2003. The status of the snow leopard in Pakistan and its conflict 
with local farmers. Oryx, 37(1): 26-33. 


Jackson R, Ahlborn G. 1988. Observations on the ecology of snow leopard 
in west Nepal. In: Freeman H. Proceedings of the 5th International Snow 
Leopard Symposium. 65-87. 

Jackson RM, Hunter DO. 1996. Snow Leopard Survey and Conservation 
Handbook. International Snow Leopard Trust, and US Geological Survey 
Biological Resources Division, Seattle, USA. 


Janečka J, Munkhtsog B, Jackson R, Naranbaatar G, Mallon D, Murphy W. 
2011. Comparison of noninvasive genetic and camera-trapping techniques 
for surveying snow leopards. Journal of Mammalogy, 92(4): 771-783. 


Janečka J, Jackson R, Yuquang Z, Diqiang L, Munkhtsog B, Buckley- 
Beason V, Murphy W. 2008. Population monitoring of snow leopards using 
noninvasive collection of scat samples: a pilot study. Animal Conservation, 
11(5): 401-411. 

Karmacharya D, Thapa, K, Shrestha R, Dhakal M, Janecka JE. 2011. 
Noninvasive genetic population survey of snow leopards (Panthera uncia) 
in Kangchenjunga conservation area, Shey Phoksundo National Park and 
surrounding buffer zones of Nepal. BMC Research Notes, 4(1): 516. 
Keehner JR. 2009. Using DNA sequence analysis of hair and feces to 
detect and identify carnivore species in National Parks of the Southwest. 
McNair Research Journal, 1: 32. 

Kimura M. 1980. A simple method for estimating evolutionary rates of base 
substitutions through comparative studies of nucleotide sequences. Journal 
of Molecular Evolution, 16(2): 111-120. 

Kohn M, Knauer F, Stoffella A, Schröder W, Pääbo S. 1995. Conservation 
genetics of the European brown bear - a study using excremental PCR of 
nuclear and mitochondrial sequences. Molecular Ecology, 4(1): 95-104. 

Liu CG, Zheng SW, Ren JR. 2003. Research foods and foods source about 
snow Leopard(Panthera uncia). Journal of 
University(Natural Science Edition), 31: 154-159. 
Long RA, MacKay P, Ray J, Zielinski W. 2008. Noninvasive Survey 
Methods for Carnivores. 2™ ed. Washington, D.C.: Island Press. 

Long RA, Donovan TM, Mackay P, Zielinski WJ, Buzas JS. 2007a. 
Effectiveness of scat detection dogs for detecting forest carnivores. Journal 
of Wildlife Management, 71: 2007-2017. 

Long RA, Donovan TM, Mackay P, Zielinski WJ, Buzas JS. 2007b. 
Comparing scat detection dogs, cameras, and hair snares for surveying 
carnivores. Journal of Wildlife Management, 71(6): 2018-2025. 


Shaanxi Normal 


Lovari S, Ventimiglia M, Minder |. 2013. Food habits of two leopard species, 
competition, climate change and upper treeline: a way to the decrease of 
an endangered species? Ethology Ecology & Evolution, 25(4): 305-318. 
Ma M, Munkhtsog B, Xu F, Mardan T, Yin SJ, Wei SD. 2005. Markings as 
indicator of snow leopard in field survey, in Xinjiang. Chinese Journal of 
Zoology, 40(4): 34-39. 

Ma M, Xu F, Munkhtsog B, Wu YQ, McCarthy T, McCarthy K. 2011. 
Monitoring of population density of snow leopard in Xinjiang. Journal of 
Ecology and Rural Environment, 27(1): 79-83. (in Chinese) 

Madden T. 2002. The BLAST Sequence Analysis Tool. /n: McEntyre J, 
Ostell J. The NCBI Handbook. Bethesda (MD): National Center for 
Biotechnology Information (US). 

McCarthy K, Fuller T, Ming M, McCarthy T, Waits L, Jumabaev K. 2008. 
Assessing estimators of snow leopard abundance. Journal of Wildlife 
Management, 72(8): 1826-1833. 


78 www.zoores.ac.cn 


Mills LS, Pilgrim KL, Schwartz MK, Mckelvey K. 2001. Identifying lynx and 
other North American felids based on MtDNA analysis. Conservation 
Genetics, 1(3): 285-288. 

Mukherjee S, Ashalakshmi CN, Home C, Ramakrishnan U. 2010. An 
evaluation of the PCR-RFLP technique to aid molecular-based monitoring 
of felids and canids in India. BMC Research Notes, 3: 159. 

Naidu A, Smythe LA, Thompson RW, Culver M. 2011. Genetic Analysis of 
Scats Reveals Minimum Number and Sex of Recently Documented 
Mountain Lions. Journal of Fish and Wildlife Management, 2(1): 106-111. 
Oli MK. 1993. A key for the identification of the hair of mammals of a snow 
leopard (Panthera uncia) habitat in Nepal. Journal of Zoology, 231(1): 71-93. 
Peres-Neto PR, Jackson DA, Somers KM. 2003. Giving meaningful 
interpretation to ordination axes: assessing loading significance in principal 
component analysis. Ecology, 84(9): 2347-2363. 

Prugh LR, Ritland CE. 2005. Molecular testing of observer identification of 
carnivore feces in the field. Wildlife Society Bulletin, 33(1): 189-194. 

R Development Core Team. 2013. R: A language and environment for statistical 
computing. R Foundation for Statistical Computing, Vienna, Austria. 

Rodgers TW, Janečka JE. 2012. Applications and techniques for non- 
invasive faecal genetics research in felid conservation. European Journal of 
Wildlife Research, 59(1): 1-16. 

Rozhnov VV, Zvychainaya EY, Kuksin AN, Poyarkov AD. 2011. Noninvasive 
molecular genetic analysis in studying the ecology of the snow leopard: 
Problems and prospects. Russian Journal of Ecology, 42(6): 439-444. 
Santini A, Lucchini V, Fabbri E, Randi E. 2007. Ageing and environmental 
factors affect PCR success in wolf (Canis lupus) excremental DNA samples. 
Molecular Ecology Notes, 7(6): 955-961. 

Shehzad W, Mccarthy TM, Pompanon F, Purevjav L, Coissac E, Riaz T, 
Taberlet P. 2012. Prey Preference of Snow Leopard (Panthera uncia) in 
South Gobi, Mongolia. PloS One, 7(2): €32104. 

Stenglein JL, De Barba M, Ausband DE, Waits LP. 2010. Impacts of 
sampling location within a faeces on DNA quality in two carnivore species. 
Molecular Ecology Resources, 10(1): 109-114. 

Sunquist ME, Sunquist F. 2002. Wild Cats of the World. Chicago, IL: 
University of Chicago Press. 

Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: 
molecular evolutionary genetics analysis version 6. 0. Molecular Biology 
and Evolution, 30(12): 2725-2729. 

Wang J, Laguardia A, Damerell PJ, Riordan P, Shi K. 2014. Dietary overlap 
of snow leopard and other carnivores in the Pamirs of Northwestern China. 
Chinese Science Bulletin, 59(25): 3162-3168. 

Wei FW, Rao G, Li M, Fang SG, Feng ZJ. 2001. Molecular scatology and its 
application--Reliability, limitation and prospect. Acta Theriologica Sinica, 
21(2): 143-152. (in Chinese) 

Xu F, Ma M, Yin SJ, Mardan T. 2005. Snow leopard survey in Tumor Nature 
Reserve, Xinjiang. Sichuan Journal of Zoology, 24(4): 608-610. (in Chinese) 
Zhang Y, He L, Duo H, Li D, Jin K. 2009. A preliminary study on the population 
genetic structure of snow leopard (Unica unica) in Qinghai Province utilizing 
fecal DNA. Acta Theriologica Sinica, 29(3): 310-315. (in Chinese) 

Zhang Z, Schwartz S, Wagner L, Miller W. 2000. A greedy algorithm for aligning 
DNA sequences. Journal of Computational Biology, 7(1-2): 203-214. 

Zuercher GL, Gipson PS, Stewart GC. 2003. Identification of carnivore 
feces by local peoples and molecular analyses. Wlidlife Society Bulletin, 
31(4): 961-970. 


